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We predict a sequence of magic-zero wavelengths for the Sr excited 5s5p 3 Po state, and provide 
a general roadmap for extracting transition matrix elements using precise frequency measurements. 
We demonstrate that such measurements can serve as a best global benchmark of the spectroscopic 
accuracy that is required for the development of high-precision predictive methods. These magic-zero 
wavelengths are also needed for state-selective atom manipulation for implementation of quantum 
logic operations. We also identify five magic wavelengths of the 5s 2 1 So~5s5p 3 Po Sr clock transition 
between 350 nm and 500 nm which can also serve as precision benchmarks. 
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The drive for increased precision in atomic, molecu¬ 
lar, and optical (AMO) measurements has led to tighter 
tests of fundamental physics m , transformational im¬ 
provements in time and frequency metrology mm, and 
suppression of decoherence in quantum information pro¬ 
cessing 0EJ. Moreover, increases in AMO precision have 
resulted in the discovery of new applications such as: lab¬ 
oratory tests of time-variation of fundamental constants 
[HHlQj, searches for topological dark matter with atomic 
clocks m and magnetometers m , probes of gravity and 
general relativity with atomic interferometry BMS3, and 
use of decoherence-free subspaces for Lorentz symmetry 
tests with entangled trapped ions [3]. 

For many of those applications, accurate knowledge of 
atomic properties has been critical for the design and 
interpretation of experiments, quantifying and reducing 
uncertainties and decoherence, and development of con¬ 
cepts for next-generation experiments and precision mea¬ 
surement techniques. Progress in development of high- 
precision theory [7511751 has yielded accurate predictions 
of many needed properties. In turn, high-precision mea¬ 
surements [19ll24j have provided experimental bench¬ 
marks for refinement and improvement of theory. 

Determination of transition matrix elements between 
excited states is a particulary difficult challenge for both 
theory and experiment. For example, the accuracy of 
theory has reached 0.2% for low-lying state transitions 
of alkali-metal atoms m and has attained 1% to a few 
percent accuracy for more complicated systems EMU- 
At present, further progress is hindered by scarcity of 
high-precision (better than 1%) benchmarks of transition 
amplitudes, which are important in many of the applica¬ 
tions cited above. Heretofore, determination of transi¬ 
tion matrix elements between excited states have been 
based primarily upon the measurements of lifetimes and 
branching ratios, which seems very difficult to push be¬ 
yond 1% accuracy, with Very few such measurements are 


available, including the most recent 3 Di lifetime mea¬ 
surement in Sr [4], where 0.5% uncertainty was achieved 
for the improved evaluation of the dynamic black-body 
shift. Moreover, such techniques face irreducible diffi¬ 
culties if the relevant branching ratios are small. This 
problem is particularly difficult for transitions involving 
excited states. 

In this work, we provide a roadmap for extracting tran¬ 
sition matrix elements using spectroscopic measurements 
of frequency. We make use of “magic-zero” frequencies 
at which the frequency-dependent polarizability, a(ui), of 
a given atomic state vanishes Emm. These magic-zero 
frequencies, which are analogues of interference phenom¬ 
ena encountered in classical coupled LC electrical cir¬ 
cuits, are ubiquitous in accessible regions of the optical 
spectrum. We demonstrate that such frequency mea¬ 
surements can serve as a best global benchmark of high- 
precision atomic theory. 

Two magic-zero wavelengths for the ground state of 
alkali-metal atoms have been recently measured, in Rb 
[29] and K [30] ■ Thus, there is established experimen¬ 
tal methodology to measure these quantities. The Rb 
measurement, which is based on finding the null point of 
diffraction of ultracold atoms by an optical lattice, was 
used to determine 5s — 6 p matrix elements as was pro¬ 
posed in [3lj . That determination attained an accuracy 
of 0.15-0.3% EH. 

In this paper, we discuss a systematic approach to the 
determination of multiple transition probabilities using 
magic-zero spectroscopy and we describe its application 
to the important test case of the 5s5p 3 Pq state of Sr. At 
present, this is one the best available benchmark systems 
for the following reasons: 

• Well-developed experimental techniques are al¬ 
ready established for Sr due to its prominence in 
atomic clock development [4j and studies of many- 
body effects in degenerate quantum gases [35, 33;. 
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FIG. 1: The dynamic polarizabilities of Sr 5s 2 and 
5s5p 3 Po states in the 350-460 nm wavelength range. 


FIG. 2: The dynamic polarizabilities of Sr 5s 2 3 So and 
5s5p 3 Po states in the 460-840 nm wavelength range. 


• The matrix element for the lowest-energy relevant 
Sr transition, 5s5p 3 Po-5s4d 3 D i, is known with 
0.23% accuracy from a recent 5s4 d 3 D i lifetime 
measurement [?]. This will simplify the extrac¬ 
tion of matrix elements from 3 Po to higher excited 
states. 

• From a theoretical prospective, Sr is one of the best 
understood atoms with more than one valence elec¬ 
tron, due to recent calculations of blackbody radi¬ 
ation (BBR) shift of the 5s 2 1 So~5s5p 3 Po atomic 
clock transition ’Tlj • 

In addition to precision measurement, there is another 
important application of the magic-zero wavelengths in 
the alkaline-earth atoms. The trapping potential in an 
optical lattice for a given atomic state is proportional to 
its dynamic polarizability a(uj). For the magic-zero wave¬ 
length, a(co) = 0, resulting in a vanishing ac Stark shift 
of that state. Thus, atoms in that state are insensitive to 
laser light of that frequency. This enables state-selective 
atom manipulation for the implementation of the quan¬ 
tum logic operations [33 ESS] - In this work, we locate 
all of magic-zero wavelengths above 350 nm for both the 
ground and 3 P 0 excited state of Sr. 

Unless stated otherwise, we use the conventional sys¬ 
tem of atomic units, a.u., in which e,m e , 4-7reo and the 
reduced Planck constant h have the numerical value 1. 
Polarizability in a.u. has the dimension of volume, and 
its numerical values presented here are expressed in units 
of cSq, where a o « 0.052918 nm is the Bohr radius. 
The atomic units for a can be converted to SI units via 
a/h [Hz/(V/m) 2 ]=2.48832xl0 _8 a [a.u.], where the con¬ 
version coefficient is 47re 0 ag//). and the Planck constant 
h is factored out. Vacuum values are reported for all 
wavelengths. 

We now demonstrate how to extract a set of matrix 


elements from a sequence of magic-zero wavelengths for 
Sr 5s5p 3 P 0 state. Consider the standard expression for 
a frequency-dependent polarizability of a state v in terms 
of the sum over all other atomic states k [57] : 


«oM 


2 ^ (k\\D\\v) 2 (E k -E v ) 

3(2J+1)^ (E k -E v ) 2 -tu 2 


(1) 


Here, J is the total angular moment of the state v and 
(k ||Z)|| v) are the reduced electric-dipole matrix elements 
which are the subject of the present work. In these equa¬ 
tions, w is assumed to be at least several line widths off 
resonance with the corresponding transitions. Linear po¬ 
larization is assumed in all calculation. 

It is self-evident from Eq. |l]) that the contribution of 
state k to the sum changes sign as lo crosses the value 
E k — E v . Magic-zero wavelengths arise due to the can¬ 
cellation of the contribution from a given resonant state 
k with the contributions from all of the other resonant 
states. Thus, there is a magic-zero wavelength between 
each pair of adjacent resonances, as can be seen in Figsjl] 
and [2] The actual location of each magic zero wave¬ 
length depends upon the distribution of electric-dipole 
matrix elements (k ||D|| v). Different matrix elements will 
be important for the determination of each magic zero- 
wavelength. Thus, measurement of a series of magic- 
zero wavelengths will enable one to extract the entire 
recommended set of the matrix elements for transitions 
to the 3 P 0 state. Furthermore, we can obtain comple¬ 
mentary information from magic wavelengths of the Sr 
5s 2 1 Sq — 5s5p 3 P 0 clock transition, at which Stark shift 
of the clock transition vanishes. With such high-precision 
benchmarks established in one system to test a new the¬ 
oretical approach, the theory can be applied to a large 
number of other systems where no experimental data are 
available. 
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Sr has two valence electrons outside a closed Kr-like 
atomic core. The main challenge in the theoretical treat¬ 
ment of systems with two or more valence electrons 
is the accurate treatment of both core-valence correla¬ 
tions and strong valence-valence correlations. We use 
the hybrid approach introduced in m that combines 
configuration interaction (Cl) and all-order linearized 
coupled-cluster methods. The core-valence (and core¬ 
core) correlations are treated by the coupled-cluster all¬ 
order method, which is used to construct the effective 
Hamiltonian. Then, this effective Hamiltonian is used 
in the configuration-interaction part of the method that 
treats the valence-valence correlations. Thus, all of the 
correlation correction to the wave functions is treated at 
the all-order level. Next, the resulting wave functions 
are used to evaluate matrix elements of various one-body 
operators, such as electric and magnetic mulitpole, mag¬ 
netic and quadrupole hyperfine, and various P-odd and 
T-odd interactions. These matrix elements are also then 
used for evaluation of polarizabilities, P-odd and T-odd 
amplitudes, and long-range interaction coefficients Cq 
and Cg. This approach is generally applicable for sys¬ 
tems with several valence electrons and has a wide range 
of applications 

To determine the magic-zero wavelengths for the Sr 
5s 2 1 S 0 ground and 5s5p 3 P 0 excited states, we need to 
calculate their frequency dependent polarizabilities for a 
wide range of frequencies. The valence part of the po¬ 
larizability is determined by solving the inhomogeneous 
equation of perturbation theory in the valence space, 
which is approximated as 

(E v - H eS )\*(v,M')} =D eSjq \* 0 (v, J,M)) (2) 

for a state v with the total angular momentum J and 
projection M |38| . While the H e ff includes the all-order 
corrections as described above, the effective dipole op¬ 
erator D e ff only includes random phase approximation 
(RPA) corrections at the present time. 

A few other small corrections can be calculated with 
second-order many-body theory. These corrections to 
dipole matrix elements contribute significantly in the 
5s5p 3 Pq BBR shift calculation of Sr [35- Their omis¬ 
sion in earlier work |39| resulted in significant difference 
of the clock dc Stark shift with experiment [23j . The cor¬ 
rections are very significant for the hyperfine constants 

US- 

The present bottleneck in the accuracy of this ap¬ 
proach can be summarized as follows. The treatment 
of corrections to the matrix elements of the one-body 
electric-dipole operator D e g , and all of the other one- 
body operators mentioned above is limited to RPA and 
second-order many-body perturbation theory. We are 
currently developing a full all-order treatment of correc¬ 
tions to one-body operators, however there are essentially 
no experimental benchmarks that we can use to estab¬ 
lish how well this approach will work. The measurements 
that we propose in this work will remedy this outstand¬ 
ing problem. We note that the method is non-specific 


TABLE I: Resonance wavelengths A and reduced dipole ma¬ 
trix elements D in Sr. Vacuum wavelength values are given 
in nm. The recommended set of the matrix elements in a.u. 
is from [34]. “Ref. [J. 


Transition Wavelength Matrix elements D 




Cl+all 

Expt. 

Cl+all 

Recomm. 

5s5p 

3 P 0 

- 5s4d 3 ZA 2642.4 

2603.1 

2.714 

2.6707(62)“ 

5s5 p 

3 P 0 

- 5s6s 3 Si 682.0 

679.3 

1.972 

1.962(10) 

5s5p 

3 Po 

- 5s5d 3 Di 484.2 

483.3 

2.458 

2.450(24) 

5s5p 

3 P 0 

- 5 p 2 3 Pi 471.5 

474.3 

2.627 

2.605(26) 

5s5p 

3 Po 

- 5s7s 3 Si 433.3 

432.8 

0.522 

0.516(8) 

5s5 p 

3 Po 

- 5s6d 3 D\ 394.3 

394.2 

1.175 

1.161(17) 

5s5p 

3 Po 

- 5s8s 3 Si 377.1 

378.2 

0.302 


5s5p 

3 Po 

- 5s7d 3 D\ 361.4 

363.0 

0.822 


5s5 p 

3 Po 

- 5s9s 3 Si 348.4 

355.4 

0.270 


5s5 p 

3 Po 

- 5s8d 3 D\ 336.7 

334.8 

0.820 



TABLE II: Magic-zero A zer o and magic A mag ic wavelengths. 
See text for the explanation of the recommended value calcu¬ 
lations. 


Cl+all-order 

Recomm. 

Cl+all-order 

Recomm. 


5s5p 

J Po 

5s 2 1 S 0 - 

5s5p 3 P 0 



355.92 


354.9 

367.0 


368.45 

358.5 

360.0 

377.8 


378.81 

376.8 

377.75 

403.35 


403.428 

390.1 

389.9 

434.35 


433.85 


497.0 

478.35 


479.126 



634.7 


632.83 



1672.9 


1666.6 




5s 2 

'So 



679.55 


689.20 




for the particular type of the one-body operator and the 
advancement in the treatment of the dipole operator will 
also provide improvement for the other operators. 

While we calculate the polarizabilities by solving the 
inhomogeneous equation |2| , which accounts for the con¬ 
tribution from all bound and continuum states, we ex¬ 
tract several dominant contributions from the low-lying 
bound states using the sum-over-states formula ([I]). To 
improve the accuracy of the calculations, we replace a few 
dominant terms in our polarizability calculations by their 
“recommended” values, which contain experimental en¬ 
ergies and recommended matrix elements from Ref. [34], 
where available. To obtain recommended values for the 
3 P 0 state, the energies of the ten lowest transitions are 
replaced by the experimental values five of the tran¬ 
sition matrix elements are replaced by the recommended 
values from Ref. [33], and 5s5p 3 Po — 5s4d 3 Di matrix 
elements is taken from [J. For the 1 Sq state, the energies 
of the four lowest transitions and the 5s 2 1 <S'o — 5s5p 1 Pi 
matrix element is replaced by the experimental values 
mi e]. Both ab initio and recommended values are 
listed in Tabled 
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TABLE III: Breakdown by transition of the contributions (in a.u.) to dynamic polarizability of 5s5p 3 Po state, at the eight 
magic-zero wavelengths indicated. The first ten rows give the contributions from the transitions indicated, and all other 
contribution are grouped together in row “Other”. The chain of dominant contributions relevant to the extraction of matrix 
elements (see text for a discussion) is highlighted in bold. 


Contribution 

1666.6 nm 

632.84 nm 

479.127 nm 433.85 nm 

403.429 nm 

378.81 nm 

368.45 nm 

355.92 nm 

5s5p 3 P 0 - 

5s4d 3 D 1 

- 188.7 

-17.1 

-9.5 

-7.8 

-6.7 

-5.9 

-5.6 

-5.2 

5s5p 3 P 0 - 

5s6s 3 Si 

45.9 

- 251.4 

-37.9 

-26.4 

-20.8 

-17.3 

-15.9 

-14.5 

5s5p 3 P 0 - 

5s5d 3 D\ 

46.3 

101.9 

-2404 

-176.0 

-97.5 

-67.6 

-58.9 

-50.3 

5s5p 3 Pa - 

5 p 2 3 Pi 

51.2 

107.5 

2361 

-241.2 

-123.2 

-82.9 

-71.7 

-60.7 

5s5p 3 P 0 - 

5s7s 3 Si 

1.8 

3.2 

9.2 

337.9 

-11.2 

-5.5 

-4.4 

-3.5 

5s5p 3 Pa - 

5s6 d 3 Di 

8.2 

12.7 

24.1 

44.6 

171.7 

-93.8 

-53.8 

-34.3 

5s5p 3 P 0 - 

5s8s 3 Si 

0.5 

0.8 

1.3 

2.1 

4.2 

147.2 

-9.5 

-3.9 

5s5p 3 Pa - 

5s7d 3 Di 

3.8 

5.4 

8.4 

12.0 

18.9 

44.0 

122.7 

-89.1 

5s5p 3 P 0 - 

5s9s 3 Si 

0.4 

0.6 

0.8 

1.2 

1.7 

3.2 

5.5 

142.7 

5s5p 3 Pa - 

5s8d 3 Di 

3.4 

4.6 

6.4 

8.1 

10.6 

15.1 

18.9 

28.6 

Other 


27.1 

32.0 

39.2 

44.8 

52.3 

64.0 

72.7 

90.6 

Total 


0.0 

0.0 

-0.3 

-0.7 

0.0 

0.4 

0.0 

0.4 

Uncertainty 

in a 

11.6 

8.7 

10.5 

6.6 

12 

67 

4.3 

2.2 

Uncertainty 

in A zero 

0.1 nm 

0.2 nm 

0.05 nm 

0.25 nm 

0.05 nm 

0.1 nm 

0.6 nm 

6 nm 


The dynamical polarizabilities of Sr 5s 2 1 Sn ground 
and 5s5p 3 Po excited states are plotted in Figs, [l] and [ 2 ] 
For clarity, we separate these into two wavelength re¬ 
gion; 350-460 nm and 460-840 nm. All resonances from 
TableHlare labeled on the top frame of the figures. For all 
of the ^Po resonances, we only write the upper transition 
states and omit 5s5p 3 Pq designations. The magic-zero 
wavelengths are determined by locating points where ei¬ 
ther 3 P 0 or 1 S'o polarizabilities vanish. We also show 
magic wavelengths, where the 3 Po and 1 S'o polarizabili¬ 
ties are identical. The 813.4 nm magic wavelength used 
for the Sr lattice clock m is shown near the right edge 
of Fig. [2] We note that our graphs show 5 other Sr magic 
wavelengths. The values of both magic-zero and magic 
wavelengths are summarized in Table [TTJ Both ab initio 
and recommended values are listed. Our value of the 
5s 2 3 So magic-zero wavelength agrees with the results of 
[35l 43] . Ref.[35] estimate of the 3 Pq magic-zero wave¬ 
length at 627 nm is 5 nm away from our 632.8(2) nm 
value. 

The contributions from the 10 lowest resonant states 
to the 5s5p 3 P 0 polarizabilities at the magic-zero wave¬ 
lengths are given in the first 10 rows of Table |III| The 
exact values of those wavelengths are given as the labels 
of the eight columns. All of the other contributions are 
grouped together in the row “Other”. This contribution 
smoothly varies with the wavelength and is substantially 
smaller than the dominant contributions for most magic- 
zero wavelengths. The sum of all contributions is zero 
within the numerical accuracy. 

Table [TTT| demonstrates the use of the magic-zero wave¬ 
length measurements to establish experimental bench¬ 
marks for the transitions noted in bold font. First 
5s5p 3 P 0 — 5s4d 3 D\ matrix element is already known 
from a precision lifetime measurement [4], so contribu¬ 
tions in the first row are known with 0.5% precision. 
The matrix element of the second transition 5s5p 3 Pq — 


5s6s 3 Si is constrained to 0.5% by the position of the 
813.4 nm magic wavelength [31]; thus the values in the 
second row are known to about 1%. 

We propose two methods for the extraction of the re¬ 
maining matrix elements. First, a global fit of all mea¬ 
sured magic-zero wavelength can be done, varying the 
dominant contributions to best match the experimental 
values of wavelength. Second, a simpler procedure can 
be used to extract the matrix elements sequentially by 
determination of dominant contributions, as follows 

1. The dominant cont ributions of the third and fourth 
transitions in Table III associated with the 5s5d 3 D\ and 
5 p 2 3 Pi states respectively, can be extracted from mea¬ 
surements of the first three magic-zero wavelengths. If all 
three are available, the uncertainty in the “Other” con¬ 
tribution can be established as well. The 479 nm wave¬ 
length is particulary useful, due to very large contribu¬ 
tions of the third and fourth transitions. The correlation 
corrections are expected to be different for the transition 
to the 5 p 2 3 P\ state in comparison to all other states, 
due to its different electron configuration (see Table II of 
Ref. [341 ). Therefore, we recommend the precision mea¬ 
surement of 479 nm wavelength as a first priority. 

2. With first four dominant contributions in hand, the 
fifth and sixth can be obtained independently from the 
434 nm and 403 nm magic-zero wavelengths since only 
one of these transitions contributes significantly to its re¬ 
spective magic-zero wavelength. This happens because of 
the larger values of the 3 D 1 matrix elements in compari¬ 
son with those of 3 Si for the higher states (see Table |I|. 

3. Then, dominant contributions for the 5s8s 3 Si, 
5s7d 3 Di, and 5s9s 3 Si transitions can be obtained if the 
last three wavelengths are known. 

The magic wavelengths listed in Table [TT] can be used as 
additional benchmarks in a similar way. 

The uncertainties in the polarizability values and re¬ 
sulting estimated uncertainties in the magic-zero wave- 
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lengths are listed in the last two rows of Table |III| The 
polarizability uncertainties are obtained by adding the es¬ 
timated uncertainties in the each of ten contributions and 
uncertainties of the “Other” term in quadrature. The rel¬ 
ative uncertainties in the polarizability contributions are 
twice the estimated relative uncertainties of the corre¬ 
sponding matrix elements listed in Table [T[ The uncer¬ 
tainties in the last four matrix elements were taken to be 
3% and the uncertainty in the “Other” contribution was 
estimated at 5%. 


In summary, we have predicted the values of nine 
magic-zero and five magic wavelengths for Sr. We demon¬ 
strated how measurements of these quantities may serve 
as a sensitive experimental benchmark for further devel¬ 
opment of higher-precision first-principles theory. 

This research was performed under the sponsorship of 
the US Department of Commerce, National Institute of 
Standards and Technology, and was supported by the Na¬ 
tional Science Foundation under Physics Frontiers Center 
Grant PHY-0822671. 
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